Electronic transport in normal-conductor/graphene/normal-conductor junctions and 
conditions for insulating behavior at a finite charge-carrier density 
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We investigate the conductance of normal-conductor/graphene/normal-conductor (NGN) junc- 
tions for arbitrary on-site potentials in the normal and graphitic parts of the system. We find that 
a ballistic NGN junction can display insulating behavior even when the charge-carrier density in 
the graphene part is finite. This efi^ect originates in the different k intervals supporting propagating 
modes in graphene and a normal conductor, and persists for moderate levels of bulk, edge or inter- 
face disorder. The ensuing conductance thresholds could be used as an electronic tool to map out 
details of the graphene band structure in absolute k-space. 

PACS numbers: 73.63.-b, 72.10.Bg, 73.63.Bd, 81.05. Uw 
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I. INTRODUCTION 



Graphene, the two-dimensional arrangement of car- 
bon atoms on a honeycomb lattice that has recently 
become available through ground-breaking fabrication 
methods^''^ possesses a wide range of unique electronic 
transport properties which originate from the conical dis- 
persion relation around the corners (K-points) of the 
hexagonal Brillouin zonei^ The low-energy theory in the 
vicinity of these points is of the form of a Dirac equa- 
tion for massless, chiral fermionsi^ The intrinsic trans- 
port properties studied on the basis of the Dirac equa- 
tion (such as the quantum Hall effect r'l^i^ii^ the minimal 
conductivity^i^a^'ii^'i^^'i^^il^ and weak localiza- 
tion corrections to the conductanc o-'^^i^^'^-'^i^^'^'^i^'^ ) there- 
fore effectively probe the graphene band structure via 
the momentum difference (5k relative to the K-points. 
On the other hand, detailed information corroborating 
the conical band structure in absolute k space has re- 
cently become accessible via angle-resolved photoemis- 
sion spectroscopyi^^iS^ 

Recent theoretical transport studies have pointed 
out that highly unconventional devices could be fab- 
ricated in patterned and gated samples of graphene, 
such as Veselago lenses^I and filters for the valley de- 
gree of freedom.^* These effects already occur in sim- 
ple, rectangular graphitic samples, so-called nanorib- 
bons, which have been studied in great detail in the 
pasti^'^i^i^i^'^i^i^i^ With few exceptions, however, 
theoretical investigations of electronic transport have 
concentrated on all-graphitic structures. In experiments, 
the ultimate electronic contacts are metallic (for illus- 
tration see, e.g., Ref. HI). Two recent works^MS 
have addressed the coupling of graphene to normal- 
conducting electrodes, in each case considering normal- 
conductor/graphene/normal-conductor (NGN) junctions 
with an armchair ribbon and zigzag interfaces, as shown 
Fig. HJa). In Ref. fiH, the graphitic part was fixed at the 
value of charge-neutral graphene, while the on-site po- 
tential in the leads was changed (the results were then 
compared to the results for a set-up in which the leads are 



also graphitic^'^). In Ref. ,4(1, the on-site potential in all 
three regions was changed simultaneously (the resulting 
finite charge-carrier density in the graphitic part greatly 
enhances the conductance of the junction). 

In this paper we systematically investigate the depen- 
dence of the electronic transport through NGN junctions 
on independent on-site potentials in the leads and in the 
graphitic region. Since the transport at finite charge- 
carrier density is anisotropic and depends on the details 
of the normal electrodes, we also consider the case of 
zigzag ribbons with armchair interfaces [Figs, [ijc, d)], 
and the case of real-space leads [Figs. [IJb, d)]. We also 
investigate how the conductance depends on bulk disor- 
der, boundary roughness and interface imperfection. 

Our results entail that a ballistic NGN junction can be 
insulating even when the charge-carrier densities in the 
leads and in the graphitic region are both finite. Con- 
ceptually, this effect can be seen as the counterpart of 
the celebrated minimal, non-vanishing conductivity ex- 
hibited by a graphene sheet at the point of nominally van- 
ishing charger-carrier densityi^iii^iiSiiiii^ii^iiiii^ii^iiiii^ 

We identify a simple mechanism for this insulating be- 
havior at finite charge-carrier density, which originates in 
the mismatch of propagating modes in the normal and 
graphitic parts of the system. In the transport across a 
ballistic interface, the transverse momentum is conserved 
(modulo reciprocal lattice vectors), and the conductance 
probes whether there are propagating modes with the 
same transverse momentum on both sides of the inter- 
face. The conductance thresholds, hence, are intimately 
related to the band structure in the normal and graphitic 
parts of the junction, which restricts the transverse mo- 
menta of propagating modes. Consequently, the conduc- 
tance thresholds could be used to deliver information on 
the band structure of graphene in absolute k-space if the 
band structure in the leads is known. Our numerical com- 
putations show that the insulating behavior persists for 
moderate levels of bulk, edge or interface disorder, and 
is only destroyed for a very rough interface. The prac- 
ticality of fabricating sufficiently clean graphene ribbons 
has been demonstrated in recent experiments^^ In other 
mesoscopic systems, the fabrication and characterization 
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of clean ballistic interfaces has reached a high level of 
sophistication."*^ 

The paper is organized as follows. Section pro- 
vides background information on the tight-binding mod- 
els used to model the NGN junctions, and on the propa- 
gating and evanescent modes in the normal and graphitic 
parts of the system. Section IIIII presents numerically re- 
sults for clean junctions and identifies regions of insulat- 
ing behavior at finite charge-carrier density. Analytical 
results are given in Section IIVI We start with an ex- 
act calculation of the conductance for the case of arm- 
chair ribbons with zigzag interfaces and lattice-matched 
leads, shown in Fig.[TJa). The results allow us to identify 
the simple mechanism for insulating behavior described 
above, which can be carried over to the other geometries 
in Fig. [TJ In Section |V] we discuss the effects of bulk, 
edge and interface disorder, as well as the effect of mode 
mixing at armchair interfaces. Conclusions are presented 
in Section IVII Appendices El (on transverse-momentum 
quantization) and|B] (on the modelling of a ballistic inter- 
face to a real-space lead) give some additional theoretical 
background on features of the tight-binding model used 
in the numerical computations. 



II. THEORETICAL BACKGROUND 

In this section we provide the theoretical background 
for the analytic calculations and numerical computations 
of the conductance of NGN junctions, which are based 
on tight-binding models and the Landauer conductance 
formula. 



A. Model Hamiltonian 

Tight-binding models of NGN junctions are shown in 
Fig. [T] The tight-binding Hamiltonian is given by 
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(1) 



where Ci is a fermionic annihilation operator acting on 
lattice site i, (ij) denotes pairs of nearest neighbors, and 
the hopping matrix elements 7^- as well as the on-site 
potential energies Vi take values as appropriate for the 
region in question. 

The graphitic region is modelled by sites on a honey- 
comb lattice with lattice constant a^'^-', hopping constant 
7*'^) and on-site potential F^'^^. The normal regions 
(r = 1,2, corresponding to the left and right leads re- 
spectively) are modelled as sites on a square lattice with 
lattice constant a^'^\ hopping constant 7'^''^ and on-site 
potential V^''\ 

In order to form an NGN junction, a graphitic region of 
length C is matched to the normal regions along graphitic 
zigzag [Figs.lUa, b)] or armchair [Figs.[Hc, d)] interfaces 
of width W. Two types of matching are considered. Fig- 
ures [Ija, c) show lattice-matched leads, where the lattice 
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FIG. 1: Tight-binding models of normal- 

conductor/graphene/normal-conductor (NGN) junctions, 
consisting of a hexagonal lattice with lattice constant 



.(G) 



coupled to square lattices with lattice constant a 



(r) 



(r-=l,2 corresponding to left and right leads respectively), 
(a) Zigzag interfaces connected to lattice-matched leads. 
Zigzag interfaces connected to real-space 



leads, a^"^^ 

lattice-matched leads, a^-^^ — a'*^' 



c) Armchair interfaces connected to 
(d) Armchair interfaces 



connected to real-space leads, a'*^' = l^'"'^''- i'^ assumed 
that the charge-carrier density in graphene and in the 
leads can be controlled independently by gates (shaded) or 
chemical doping, which shifts the on-site potentials 1/'*^' (in 
the graphitic region), V"(i) (in the left lead) and T/'^' (in the 
right lead). 
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constant is related to the lattice constant in graphene 
by a'^^'^ — -\/3a^'^^ and a^^^ = a^^\ respectively. Figures 
[UJb, d) show real-space leads, approximated by a finer 
lattice with a reduced lattice constant. The hopping con- 
stants across the right and left interface are denoted by 
^(iG) gj^f^ 7*-^'"'-', respectively. 



B. Landauer conductance formula 

For small bias voltages the phase-coherent conductance 
of a mesoscopic structure is given by the Landauer for- 
mula 

g = {2e'^ lh)iYth, (2) 

where t is the transmission matrix with elements tnm- 
The mode index n refers to incoming propagating modes 
in the left lead, while the mode index m refers to outgo- 
ing propagating modes in the right lead. The interface 
couples these modes to the propagating and evanescent 
modes in the graphitic scattering region. The remain- 
der of this background section compares the properties 
of these modes in the normal and graphitic parts of the 
NGN junctions. 



C. Dispersion relations 



For the zigzag orientation of the interface, the ampli- 
tudes on the A and B sites are related via V'b(0,0) = 
'7]j|-'0a(O, 0), where the function 

/(fc.,fcj,) =7^'^^[l + 2e'''==='^"''/'cos(V3feya(«V2)] (6) 

also delivers the graphene dispersion relation via E — 
y{G) „ 77I/I. The index 77 = ±1 distinguishes the two 
branches of the dispersion relation. For the armchair in- 
terface, the graphene lattice is rotated by 90°. The am- 
plitudes 7/15(0,0) = ryjjj-VjyilO, 0) and dispersion relation 

E — V'*^) — rj\ f \ are then determined by the hmction 

The quantization of the transverse momentum ky in a 
wire geometry is discussed in Appendix \K\ 

D. Mode characterization 

Whether a mode with given transverse momentum is 
propagating or evanescent is determined by the disper- 
sion relation of the region in question. For a given trans- 
verse momentum, the dispersion relation delivers the lon- 
gitudinal wavenumber kx as a function of energy and on- 
site potential. A mode is propagating when kx is real 
and evanescent when kx is complex. 



The properties of the modes in the normal and 
graphitic regions follow from the dispersion relations, 
which relate the wave number k = {kx, ky) to the energy 
i?(k) of Bloch waves and reflect the symmetry properties 
of the underlying lattice structure. 

The unit cell of the hexagonal lattice contains two in- 
equivalent sites A and B with different orientation of the 
bonds. The unit cell of the square lattice contains only 
a single site, which we denote by S. We use the symbols 
ipAix, y), i^six, y), and tpsix, y) to denote the amplitudes 
of the wavefunction on each site, where cc, y are the co- 
ordinates of the center of the unit cell. We assume that 
one unit cell is centered at the origin x — y = Q. 

The square lattice supports Bloch waves 



V's(x,y)= 7/5(0,0)6^'===-+^'=='^ 



(3) 
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with dispersion relation 

E = Vt*^' - 27('')[cos(a('')/cj,) + cos{a^''^ kx)]. (4) 

In the continmmr limit a^"^^ at fixed ^^'^Ha^'^^)'^ = 7?—, 
one recovers the parabolic dispersion relation 

E = V^^^-ij^^^ + ^ikl + kl). (5) 

The hexagonal lattice supports Bloch waves of the form 
{',PA{x,y),^PB{x,y)) = (7Aa(0, 0), 7/5(0, 0))e^'=^-+^'=»^ 



FIG. 2: (Color) Dependence of the Fermi lines on the on-site 
potential 1/'*^' or l^'*"' in the Brillouin zones of a hexagonal 
lattice with hopping constant 7'*^' (middle panels), a lattice- 
matched square lattice with hopping constant 7'""' = 7'*^' (left 
panels) and a real space lattice with lattice reduction factor 
q = 5 and hopping constant 7''^-' = g^7''^-' (right panels). In 
the upper panels, the lattice constants are related as in the 
NGN junctions with zigzag interfaces, shown in Fig. [T][a, b). 
In the lower panels, the lattice constants are related as in the 
NGN junctions with armchair interfaces, shown in Fig. [ijc, 
d). The arrows indicate how the Brillouin zones have to be 
matched to form the four different NGN junctions in Fig. [T] 
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Propagating modes at the Fermi energy Ep = can 
be identified fi'om the condition that the fine of constant 
ky crosses one of the Fermi hnes, which depend on the 
on-site potential as shown in Fig. (5] This dehvcrs the fol- 
lowing conditions for propagating modes on the various 
types of lattice: 



27(G) 



< 



cos 




< 



G)| 



27(0) 



(hexagonal lattice with zigzag interfaces), 



(8a) 



< 



cos I -kyU^'^^ 



for <7(«) 



< 



for > V5j 



(G) 



(hexagonal lattice with armchair interfaces), (8b) 



cos 



27('-) 
(square lattice) 



< 1 



For each type of lattice, these conditions deliver the range 
of transverse momentum in which the modes are propa- 
gating, while in the complementary range the modes are 
evanescent. The border between these ranges defines the 
threshold values of transverse momentum at which the 
modes change their character. In Section IIVI we trans- 
late these thresholds into conductance thresholds for the 
NGN junctions. 
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FIG. 3: (Color) Dependence of the conductance g (in units of 
the quantum of conductance go ~ 2e^ /h) on the on-site poten- 
tials yt'^) in graphene and V'-'') in the leads. The conductance 
is color-coded as indicated in the scale bar, and the white con- 
tour denotes the conductance thresholds where g = The 
panels (a-d) correspond to the NGN junction shown in Fig. 
[ija-d). The width and length of the graphene region were 
chosen to be W = 175 a*''-' and C = 150 a^'^K Further param- 
eters for the plots are described in Section [III Al 



III. NUMERICAL RESULTS FOR CLEAN NGN 
JUNCTIONS 

A. Method and parameters 

In order to obtain an immediate insight into the gate- 
voltage dependence of the conductance of the NGN junc- 
tions of Fig. [T] we first present numerical results for clean 
systems. 

The transmission coefficients tnm are computed em- 
ploying an efficient decimation scheme4^ In this scheme, 
the uncoupled Hamiltonians of the leads and graphene 
are reduced to effective Hamiltonians (self-energies) at 
the interfaces. For square lattice leads the self-energies 
are known analytically^^ In the graphitic region, the 
renormalisation procedure is performed iteratively site- 
by-site employing Gauss- Jordan elimination. The Dyson 
equation is then used to determine the surface Green 
function of the coupled system. Finally, the transmis- 
sion coefficients follow from the Fisher-Lee relation4^ 

In our numerical computations, the graphitic region 
has width W = 175 a^*^) and length C = 150 a^"^). In 
the case of lattice-matched leads, panels (a) and (c) in 
Fig. [U the hopping constants in the leads and across the 
interfaces are all taken to be identical to the hopping 



constant in graphene, 7'''^ = 7(''G) = 7(G) ^ correspond- 
ing to a ballistic interface (without a tunnel barrier). In 
order to model real-space leads, panels (b) and (d) in 
Fig. [U the lattice constant is reduced by a factor q = 5 
as compared to the lattice-matched leads. The hopping 
constant 7^'"^ = q^7'-*^^ in the leads is chosen to preserve 

the effective mass m — / a^^'^ ] in the parabolic 
region of the dispersion relation at the bottom of the 
band. The value 7('"G) = 7.8617^'^' for the interface hop- 
ping term is again chosen to model a ballistic interface 
(for a derivation of this value see Appendix |BJ . 



B. Results 

Figures[3I^a-d) show the gate voltage dependence of the 
conductance of clean NGN junctions, where each panel 
corresponds to one of the configurations in Fig. [Ija-d). 
The results are presented in a color scale where red cor- 
responds to a large conductance, while blue corresponds 
to a low conductance. The unit of conductance is the 
conductance quantum 170 — {2e^ /h). Regions of low 
and high conductance are separated by a white con- 
tour at 5 = igo- The on-site potential T^(G) jjj ^j^g 
graphitic part is varied independently of the on-site po- 
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tential V^^^ = V^^) in the leads. These energies are mea- 
sured in units of the hopping constant in graphene, 7^'^-'. 
For lattice-matched leads, the on-site potentials are var- 
ied over the complete bandwidth of the dispersion rela- 
tion in the graphitic region (—87*^'^^ < V^'^'' < 87^'-^^) 
and in the leads {~4:^^^^'> < V^'''> < 47*'^)). For real- 
space leads the on-site potential is restricted to the range 
_4^(G) _^ 4g2^(G) < yC'-) < 4^(G) _^ 4g2^(G) in the 

parabolic region at the bottom of the square-lattice dis- 
persion relation, where g = 5 is the lattice-constant re- 
duction factor introduced above. 

The results in Fig. [3] show a highly systematic de- 
pendence of the conductance on the on-site poten- 
tials. The conductance is small for V"'*^' = 0, 
the region of minimal conductivity at zero charge- 
carrier density discussed in previous transport studies 
of graphene.^'i'^iia'iiii^'i^ii^^'i^i^'i^ However, we also 
find regions of very small conductance where the charge- 
carrier density in graphene and in the leads is finite. 
Conditions for this insulating behavior are determined 
in Section [IVl The general raggedness of the contours of 
constant conductance is a common feature in graphene 
transport studies, and can be associated to Fabry-Perot 
resonances i36i37,4o 

For lattice-matched leads with zigzag interfaces 
[Fig. inia)], we observe an approximate mirror symme- 
try of the conductance for (F^.F^^)) (V"W, 
This symmetry is most pronounced for the regions of low 
conductance, delimited by the white threshold contour. 
On the other hand, the maximal conductance for posi- 
tive T/^*^) is much larger than for negative V^^\ These 
maxima are found at V^^^^ « and F^'-^) « ±7^''^. 

For lattice-matched leads with armchair interfaces 
[Fig. M^c)], the mirror symmetry {V^''\V^^^) 
(yir) ^ —Y(.G)"j is only observed for the region of low 
conductance with y'-'"' « —47^'-'', close to the top 
of the square-lattice dispersion relation. The region 
of high conductance obeys an approximate symmetry 

when both on-site potentials are inverted, {V'^''\V'^'^^) 
(_yM,_^(G)). 

For real-space leads [Fig. [3]^b, d)], the general features 
of the conductance are inherited from the behavior for 
the lattice-matched leads in the parabolic region at the 
bottom of the square-lattice dispersion relation. In this 
region, the conductance in general increases for increas- 
ing charge-carrier density in the leads (corresponding to 
smaller values of y'-'"'), and large values of the conduc- 
tance are predominantly found for positive V^'^\ Despite 
being more ragged, the threshold contours have a similar 
general trend as in the lattice-matched case. For zigzag 
interfaces, the region of insulating behavior obeys the ap- 
proximate mirror symmetry, while this is not the case for 
armchair interfaces. 



IV. ANALYTICAL RESULTS 

Most of the conductance thresholds observed in the 
numerical computations. Fig. [31 can be explained via a 
simple mechanism, based on the mismatch of propagat- 
ing modes on both sides of an NG interface. We start our 
considerations with the exact calculation of the conduc- 
tance of ballistic NGN junctions with zigzag interfaces 
and latticed-matched leads [Fig. [H^a)]. The calculation 
shows that in this case, the insulating regions correspond 
to conditions where the propagating modes in the normal 
leads only couple to evanescent modes in the graphitic 
scattering region. This observation is then carried over 
as a criterion to calculate conductance thresholds for the 
other three types of NGN junctions. 



A. Conductance for NGN junctions with zigzag 
interfaces and lattice-matched leads 

For ballistic NGN junctions with zigzag interfaces and 
latticed-matched leads [see Fig. [Ija)] the conductance 
can be calculated analytically via a wave matching pro- 
cedure. The calculation succeeds because for the zigzag 
configuration the hard-wall boundary conditions in the 
N and G parts select the same transverse wave compo- 
nents, Eqs. (|AllA2p . Hence no mode mixing occurs at 
a clean zigzag interface. The transmission matrix be- 
comes diagonal, and the wave matching for each fixed 
transverse-mode profile n reduces to a one-dimensional 
problem. 

A derivation of the matching conditions for the 
present geometry has been given in Refs. [solliol . 
The wavefunction in the square leads, ipsix^y) = 
'tps{0,0)e^''''^'^^''^y , has to be matched with the wave- 
function in the graphitic region, {ipA{x,y),ipB{x,y)) — 
(i/;A(0,0),?/'B(0,0))e^*=-^+^'=«2' at the interfaces of the 
NGN junction (located at a; = and x = L), subject 
to the boundary conditions 

7(^^^(0,0) = 7(i^)7/.s(0,0), 

7(^'=Vb(0,0) = 7''Vs(a('\0), 

7(2<=Va(L,0) = 7''Vs(i-a('\0), 

7(2)v;b(L,0) = 7(2G)^5(L,0). (9) 

In Refs. [39II40I these equations have been solved for the 
cases y^*^) — 0, V*^^^ — V^^'' of charge-neutral graphene 
and y*-^-* = y^^-* = y^*-^) for uniformly gated junctions, 
respectively. For the general case of independent on-site 
potentials and coupling constants we find 
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-4zCsin (^|a(«)fcf ) j sin (V3a(«)fc(i)) sin (VSaC^)^- 

1 y(G)p2 y(G) 



riA(i)A(2) 7(G)A(1) 7(G)r2A(2) 



sm 



+ 2CTi sin 



£ _ ^a(G) ) fc(G) 



2C 



riA(i)A(2) " 

where the n dependence of the above expression is imphcit in both and ky , and 



C+ta(o)]ki^^ 



C = cos(^a(«'fc, 



^(1G)^(2G) 



^^(1)^(2)^(G) 



7(20) 



(10) 



(11) 



The conductance of the junction follows from the Lan- 
dauer formula In the limit of V^'^'' = and 

y^^' — y'^' one recovers the result of Ref. [s^, while for 
yii) ^ y(2) ^ y(G) ^j^g j-gg^j^ of Ref-Ois obtained. We 

confirmed that the values of conductance obtained from 
Eq. ([To]) for general combinations of the on-site potentials 
are in numerical agreement with the results obtained in 
the previous section [Fig. [3{a)] . 



B. Conductance thresholds in large ballistic NGN 
junctions due to mode mismatch 

For a long graphitic region, with C ^ a^^\ a transmis- 
sion coefficient ^U)) tends to zero when the propagating 
mode in the lead couples to an evanescent mode in the 

(G) 

graphitic region (i.e., if Imfci ^0). The conductance 
of the junction is small if this is the case for all trans- 
mission coefficients. A small conductance, hence, does 
not necessitate that all the modes in the graphitic region 
are evanescent ~ it suffices that the propagating modes 
in the graphitic region do not couple to the propagat- 
ing modes in the leads. Consequently, the conductance 
can be small even when the charge-carrier density both 
in the graphitic region as well as in the normal leads is 
finite. We now apply this mode-mismatch mechanism to 
calculate conduction thresholds of long and wide ballistic 
NGN junctions with clean interfaces, covering all of the 
geometries shown in Fig. [1] 

The requirement of a wide junction arises from the fact 
that the mechanism described above relies on the con- 
servation of transverse momentum (modulo reciprocal- 
lattice vectors). For a clean zigzag- interface, this con- 
servation law is exact. As discussed in Appendix \^ 
for an armchair interface, the quantized transverse mo- 
menta in the graphitic part, Eq. (jA3[) . differ from the 
quantized transverse momenta in the normal leads, Eq. 
(IA2p . The resulting mode-mixing for interfaces of finite 
width is discussed in Section IV Bl Moreover, details of 
the transverse-momentum quantization in graphene rib- 



bons depend on the chemistry of the edgesj 1 ' In wide 
junctions the modes in N and G are tightly spaced and 
can be assumed to be quasi-continuous. In the limit 
W CX3, the transverse wave number ky is then con- 
served exactly, also for armchair interfaces. 

The requirement of a long, ballistic NGN junction 
is needed so that we can use the assumption that the 
evanescent modes only give negligible contributions to 
the conductance, as explained in more detail below. (The 




-4 4 92 100 



(c) v<V°' V<V°' 

FIG. 4: (Color online) For each of the NGN junction shown 
in Figs. [TJa-d) , the white regions denote conditions for insu- 
lating behavior due to the mismatch of propagating modes. 
The boundary of the regions are the conductance thresholds 
derived in Eq. (|12p (for panel a), Eq. (|13p (for panel b), and 
Eq. (|14|l (for panel c). No regions of mismatching modes oc- 
cur in panel (d). Comparison to Figs. (J^b, d) shows that the 
conductance can be small even beyond the mode-mismatch 
mechanism. 
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effects of disorder are discussed in Section IVl) 

Under these conditions, good conduction requires val- 
ues of the on-site potentials V^'^^ and F'""' at which 
one finds transverse momenta, possibly differing by 
reciprocal-lattice vectors, for which the associated modes 
in N and G are both propagating. The conductance is al- 
ways small when this criterion is not fulfilled. The thresh- 
old values of the on-site potentials separating the regions 
of matching and mismatching propagating modes can be 
derived from Eqs. ([5]). 

For the arrangement in Fig. [1] (a), the conductance 
threshold due to mode mismatch has two branches given 
by 



y('") 



7 



(r) 



\ 7(G) 



1 



(12a) 
(12b) 



For a finer discretization of the square lattice as in 
Fig. [1] (b), the first branch remains within the region 
of the parabolic dispersion at the bottom of the band, 
while the second branch shifts to F^'-* —00. In the 
continuum limit ([5]) of the dispersion relation, conduc- 
tance thresholds due to mode mismatch exist only for 

■y/2m(47('-)-'V('-)) 



-\/3a^'^-' < 27r/3, and then are given by 



IfC^^I _ 2cos ' ^2^(47^-1^^) , 



1 



(G) 



(G) 



1. (13) 



The survival of conductance thresholds in the contin- 
uum limit can be best understood by considering the 
modes with fcy « 0, which propagate in the normal 
leads for on-site potentials close to the bottom of the 
parabolic dispersion relation. For zigzag interfaces, these 
modes propagate in the graphitic region only in the re- 
gion 

|y(G)| > ^(G) [ggg Fig. 12]. 
For the arrangement in Fig. [T] (c) , care has to be taken 
for the fact that the periods of the Brillouin zones of the 
leads and the graphene part differ by a factor of 3/2 in 
the ky direction (see the alignment of the Brillouin zones 
in Fig. [2]). The graphitic region hence mediates the cou- 
pling of lead modes with different transverse momentum. 
Propagating modes always match up for v = > — 3, 
while in the region v < —3 there are three branches of 
conductance thresholds due to mode mismatch. Two of 
these branches are bounded by the condition 



v(G) 



9 3 1 
'5 + 4\/l + -u + -u2 -f -w3, (14a) 



The third branch is given by the condition 



v(G) 



9 3 . 

-—V 

4 2 



(14b) 



For a finer discretization of the square lattice as in 
Fig. [l](d), all these branches move to V^'^'' —00. This 



results in the absence of conductance thresholds due to 
mode mismatch in the continuum limit. This can be 
understood from the observation that in graphene with 
armchair interfaces, one can find propagating modes with 
fcj^ « for all values of the on-site potential F'-*^^ [see 
again Fig. [2]. 

Figure [4] shows the conductance thresholds due to 
mode mismatch in the l/''^-' - F^'"-' plane. Each panel 
corresponds to one of the various types of NGN junctions 
shown in Fig. [T] Comparison with Fig. [3] shows that for 
zigzag interfaces [panels (a) and (b)], the mode-mismatch 
mechanism explains all conductance thresholds. For arm- 
chair interfaces, the mode-mismatch mechanism explains 
the conductance thresholds in the left part of panel (c), 
corresponding to energies close to the top of the band 
in the leads. The numerical results in Fig. [3Jc, d) ex- 
hibit additional thresholds in the lower-right corner of 
the plane, corresponding to the bottom of the 

conduction band in the leads and the top of the conduc- 
tion band in the graphitic part of the system. Here the 
propagating modes on both sides of the interface differ 
drastically in their longitudinal wavenumber (and hence 
in their self-energy) , which also inhibits their coupling.— 
Consequently, the conductance of a ballistic NGN junc- 
tion can be small even for conditions not described by 
the simple mode-mismatch mechanism. 



C. Sharpness of thresholds 

Above we have ignored the contribution of evanescent 
modes in graphene. These modes become important for 
a finite system size L, and determine the sharpness of the 
conductance thresholds. 

The role of the evanescent modes is best understood 
by considering the most slowly decaying modes in the 
graphitic region, and in particular by investigating which 
of these modes still couple to propagating modes in the 
leads when one enters the insulating regime. It comes in 
handy that the most slowly decaying modes have trans- 
verse wave numbers just at the threshold to where such 
modes become propagating, which is determined by Eq. 
([5]). Inside the insulating region, not only all the prop- 
agating modes in graphene, but also the adjacent slowly 
decaying evanescent modes couple to evanescent modes 
in the leads and hence do not contribute to the trans- 
port. The remaining graphitic evanescent modes which 
do couple to the propagating modes in the leads all have 

(G) 

a finite decay constant Im fci > k, and their total con- 
tribution to the conductance is suppressed exponentially 
with exp(— kX). The sharpness of the thresholds hence 
increases exponentially with the system size. 

The decay constant k approaches zero as one ap- 
proaches the conductance thresholds. Let us assume 
that this is induced by changing the on-site potential 
]/('■) ^ ^(r),thresh g^gj y (G) ^ ^^^^.^ depending on the 
geometry yC''). thresh determined by Eq. dH]), (US]), or 
(1141) . For a vanishing charge-carrier density in graphene. 
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(c) y(')/y(G) y(')/y(G) (d) 

FIG. 5: (Color) Same as Fig. |3ja), but in the presence of 
bulk disorder with u = 0.2 (panel a), surface roughness with 
/ = 0.3 (panel b), and interface disorder with F — 0.17''^^' 
and r = 7''"'^' (panels c and d, respectively). 

= 0, the linear dispersion relation close to the Dirac 
point then entails that k oc jF^'') — y(»')ithrcsh| increases 
linearly with the distance to the threshold, while for a fi- 
nite V'-'^^ it increases faster, as k cx _ yC''), thresh 1 1/2^ 
Hence, the conductance thresholds are sharper at a finite 
charge density. 

It is insightful to contrast the exponential suppression 
of the conductance carried by evanescent modes in the in- 
sulating region with their contribution inside the conduc- 
tive region. In this case the most slowly decaying modes 
do couple to propagating modes in the leads. For a van- 
ishing charge-carrier density in graphene, I/^'^' = 0, the 
conductance carried by the evanescent modes adds up to 
a contribution oc W/L, which is constant and finite at a 
fixed aspect ratio even when the system is very largeJ^i^S 
For finite V^'^\ on the other hand, their contribution is 
proportional to M^/L"' for L > a'^^^j'^'^y\V^'^'> \ and hence 
decays algebraically for increasing system size.'*° Both ex- 
pressions require that the most slowly decaying evanes- 
cent modes in the graphitic part couple to the propa- 
gating modes in the leads, which is not the case in the 
insulating regime. 

V. MODE MIXING 

The derivation of conductance thresholds from the con- 
dition of mismatching propagating modes in Section IIVI 
relied on the conservation of transverse momentum. In 
this section we explore how violations of this assumption 



modify the threshold conditions. 

A. Mode mixing by disorder 

In order to explore the effects of mode mixing by disor- 
der, we implement three different scattering mechanisms: 
short-ranged bulk disorder and surface roughness in the 
graphitic region, as well as imperfections at the NG in- 
terfaces (long-ranged bulk disorder does not provide ef- 
ficient mode- mixing). Bulk disorder is modelled via a 
random on-site potential Vi — V^'^^ + Ui where the Ui are 
independently and identically distributed (i.i.d.) random 
numbers drawn with uniform probability from an inter- 
val [— m/2,u/2]. For a rough edge we randomly elimi- 
nate a fraction / of the graphene sites within a distance 
of 20^*^^ from the boundaries of the system. An imper- 
fect interface is modelled via random hopping elements 
lij — 1^^'^'^ ij for the links crossing the interface, where 
the Tij are i.i.d. random numbers drawn with uniform 
probability from an interval [— r/2,r/2]. 

Figure [5] presents the results for an NGN junction with 
zigzag interface and lattice-matched leads [the geometry 
of Fig. [Ija); the results for the other geometries are sim- 
ilar]. Panel (a) shows the conductance for bulk disorder 
of strength u ~ 0.27'^'^^ Panel (b) shows the conduc- 
tance for surface roughness with / = 0.3, the value for 
which we found the strongest effect on the conductance. 
In both cases, the maximal conductance is reduced to 
about 2/3 of the value found in the clean case. This is 
comparable to what is found in other transport studies 
at similar levels of disorder . In contrast, the 
threshold contours delimiting the region of low conduc- 
tance are only weakly affected by the disorder. 

Figure EJc) shows the conductance for an imperfect in- 
terface r = 0.l7'-'''^-'. This moderate level of imperfection 
has only a minimal effect on the regions of high and low 
conductance. A noticeable change of these regions is only 
induced for a rough interface, shown in Fig. El^d), where 
the fluctuations F = 7''"'^' are set equal to the average 
interface hopping element. At this level of imperfection, 
the regions of low conductance cover a smaller part in 
the V^^'^ - plane. This has to be attributed to the 
diffractive effects of a rough interface. It is interesting to 
observe that the conductance thresholds are most robust 
around V^'^^ — 0; especially, the conductance threshold 
in the region V'^'^^ < is almost unchanged even though 
the interface is very rough. 



B. Mode mixing at clean armchair interfaces 

As discussed in Appendix \K\ the quantized transverse 
momenta in the graphitic part of an NGN junction with 
armchair interfaces (and zigzag surface), Eq. (|A3|) . dif- 
fer from the quantized transverse momenta in the nor- 
mal leads, Eq. (|A2p . This results in a finite amount 
of mixing even for a clean interface, which is automat- 
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to the charge-neutrality point, where the charge-carrier 
density nominally vanishes .^ii^iiaiiiii^ii^iiiii^i^iili^ 

We identified a simple mechanism for the conductance 
thresholds at finite charge-carrier density, namely the de- 
coupling of the propagating modes in the different parts 
of the system due to the mismatch of their transverse 
momenta. Since these momenta are determined by the 
dispersion relation, the conductance thresholds could in 
principle be used to obtain information about the band 
structure of graphene if the dispersion relation in the 
leads is known. Our numerical computations show that 
such an analysis would be robust against the effects of 
bulk and surface disorder in the graphitic region, and 
would also tolerate a moderate amount of imperfection 
of the interfaces. 



FIG. 6: Density plot of the transmission coefficients |tnm| 
for an NGN junction with armchair interfaces and lattice- 
matched leads [the geometry of Fig. [ijc)], with 1/'*^' = 7'"^^ 
V^^^ = (other parameters as in Fig. |3jc)). 



ically accounted for in the numerical results of Section 
ITTl Figure [6] shows a density plot of the modulus \tnm\ 
of the transmission amplitudes for the NGN junction in 
Fig. [ijc), with parameters as for the computations in 
Fig. EJc). The on-site potentials are set to the values 
yiG) = 7^'^), = 0, where all modes are propagating, 
so that the mode mixing can be seen in transmission. The 
figure shows that the transmission matrix is sparse. Each 
mode mixes with a small number of modes with similar 
transverse momentum. The additional branches origi- 
nate from the different periodicity of the Brillouin zones 
in the leads and the graphitic part, already mentioned 
in Section [TVl At each interface, the transverse momen- 
tum is only conserved modulo reciprocal lattice vectors. 
The periodicity of the Brillouin zone of the leads and the 
graphitic region in ky direction differs by a factor of 3/2. 
For the conditions of Fig. [51 this mediates the coupling 
into two additional branches of transverse momenta in 
the leads. 



VI. CONCLUSIONS 

In this work we systematically investigated the gate- 
voltage dependence of the conductance of four vari- 
ants of normal-conductor/graphene/normal-conductor 
(NGN) junctions, consisting of a graphene strip which 
is coupled in different ways to normal leads of identical 
width (c.f. Fig. [T]). Starting from exact numerical com- 
putations we identified conditions of insulating behavior 
in clean junctions, which can be encountered even when 
the charge-carrier density in the central graphitic region 
and in the normal-conducting leads is finite. Concep- 
tually, this effect can be seen as the counterpart to the 
celebrated minimal, finite conductivity of graphene close 
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APPENDIX A: TRANSVERSE-MOMENTUM 
QUANTIZATION 

In a wire geometry, the boundary conditions select 
a discrete set of W transverse wave numbers ky for a 
given propagation or decay direction, which we enumer- 
ate by a mode index n = 1,2,3,..., W. The details of 
the transverse-momentum quantization of graphene rib- 
bons depend on the chemistry of the edgesj^i^i^ In this 
paper, we are mostly concerned with wide graphitic re- 
gions, where the transverse momentum becomes quasi- 
continuous. When we, in the following, give expressions 
for W in the tight-binding model used in the numerical 
simulations, it should be noted that the dimension W 
refers to the width of the interface, which is identical to 
the width of the square-lattice leads but differs from the 
width of the graphitic region (which is wider by VSa^*^^ 
for zigzag interfaces and by a^*^-* for armchair interfaces). 

On the square lattice, W =1 + W/a^'') is the number 
of sites in the cross-section of the wire, and the set of 
quantized transverse wave numbers is given by 



k 



(Al) 



On the hexagonal lattice with zigzag interfaces (and 
hence armchair boundaries), W ^ 1 +W/^/3a'-'^\ and 
the set of quantized transverse wave numbers is given by 



ky = -. A2 

Both sets of quantized transverse wave numbers become 
identical when the lattices are matched commensurably 
as shown in Fig. [IJa), where a*^''^ = ^/3a^'^\ 
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More complicated is the case of the hexagonal lattice 
with armchair interface (which has zigzag edges), shown 
in Fig.[Hb). At the upper (lower) edge, hard- wall bound- 
ary conditions translate into a vanishing amplitude on 
the A (B) sites in the first unit cell beyond the wire 



boundary. Since the amplitudes on these sites carry a rel- 
ative phase which depends on the propagation direction, 
the quantized transverse wave numbers are determined 
by a transcendental equation, 



^^^(0)2/^(02 _ sin2(3fcj,a(G')/2) + i sin(3fcj,a(«) /2) 
^YiG)-^ l^(G)'^ - sin2(3fcj,a(G)/2) - i sin(3fcj,a(G) /2) 



(A3) 



r 



left lead 



right lead 
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FIG. 7: Interface between two commensurably matched 
square-lattice leads of different lattice constant (lattice- 
constant reduction factor q = h\ 



For < 7*^*^^ and close to the K points, this equa- 

tion reduces to the condition derived from the Dirac 
equation given in Ref. [s^ In general, Eq. (IA3P has 
W = (4/3) 2>V/3a('^) independent solutions. For 
I ^ 7'*^), this includes a number of edge states with 
Imfcj, 7^ 0. For wide interfaces (VF ^ 1), the real- valued 
transverse wave numbers are almost uniformly spaced, 
but do not coincide with the transverse wave numbers 
of the lattice-matched square lattice with a'^''^ — 0^*^^ 
shown in Fig. [T]Jb) . 



APPENDIX B: INTERFACE HOPPING 
CONSTANT FOR A TRANSPARENT 
INTERFACE WITH A REAL-SPACE LATTICE 

In this appendix we describe how to model a trans- 
parent interface between a tight-binding lattice and a 
real-space lattice, as shown in Fig. [Hb,d). A simpler 
variant of such an interface, well suited for analytical cal- 
culations, is the interface of two commensurably matched 
square-lattice leads with different lattice constant, shown 
in Fig. [T] The lattice constants in the left and right lead 
are related by an integer-valued reduction constant q so 



that a(^) 



lq\ the continuum limit of the right lead 



representative for couplings of real space leads to other 
tight-binding lattices, since only the sites adjacent to the 
interface enter the subsequent considerations. 

The hopping constants in the left and right lead are 



denoted by 7^-^^ and 



(R) 



respectively. In our nu- 



merical computations we further made the assumption 
that these constants are related by the requirement of 
an identical effective mass to = f? I\l'^'^^^{a^^'^')'^\ = 
/i2/[27(^)(a'^')2]. The hopping constants in the two 



leads are then related by 7''''^'' 

Our goal is to determine the inter-lead hopping con- 
stant, 7^^^), so that the lead is transparent for energies 
in the parabolic region of the dispersion relation at the 
bottom of the bands of both leads. 

We choose a coordinate system where the right lead 
starts at a; = and the origin accommodates one of the 
lattice sites that is linked to the left lead. Now consider a 
particle arriving from the left lead at a fixed energy E and 
transverse wave number ky^ which are both conserved 
under reflection from the interface. The wave function 

V;(^)(x, y) = Ae'^'i'^ ^"^""^^y -I- Be-*'=-'''^+^'=''^ (Bl) 

then describes the superposition of the incoming and re- 
flected particle, where the longitudinal wave number fci^^ 
is fixed by the dispersion relation ([3]). 

Upon crossing the interface, the transverse momentum 
is conserved modulo reciprocal lattice vectors. Since the 
Brillouin zone in the right lead is larger by a factor 5, 
one couples into q inequivalent modes with wavenumber 
^y,v — + 27rp/a^^\ p — 0, l,...,g — 1. For energies 
close to the bottom of the bands in both leads, only the 
mode with p = is propagating while the others are 
evanescent. We denote the longitudinal wavenumber of 
the propagating mode by fci , and the decay constant 
of the evanescent modes by Kp, p — 1,2,. — 1. The 
wave function in the right lead is hence given by 



9-1 

p=l 



z+2'Kipy / a 



is approached for q ^ 00. This simple arrangement is 



(B2) 

The relative amplitudes Cp of the evanescent modes 
are fixed by the boundary condition on those sites at 
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a; = which have no link to left lead. Since the 
wavefunction (|B2p would fulfil the Schrodinger equa- 
tion when the right lead would be continued beyond the 
interface, this boundary condition can be formally ex- 
pressed as 'ilj^-^\—a^-^\y) = 0, where y = a^^'^p', with 
p' = 1, 2, . . . , (7 — 1, is the transverse coordinate of in- 
equivalent disconnected sites. These boundary condi- 
tions are fulfilled when the amplitudes take the value 

Cp — c 

The remaining amplitudes A, B, and C are now ob- 
tained from the boundary conditions of the sites on both 
leads which are linked to the other lead. Using again the 
fact that the wavefunctions (jBip and (jB2p both fulfil the 
Schrodinger equation when the leads would be extended 
across the interface, these conditions can be written as 

^(flX)^(i?,)(O^Q) ^^(L)^(L)(Q^O)^ (B3) 
^W^(fl)(_a(^),0) =7(^^)V^'^n-«^^\0). (B4) 
The reflection amplitude r = B/A then follows as 

where 

g = 1 + ^ Cp (B6) 



At the bottom of the band, we can further assume 
ky,k^'^^ « 0, and the decay constants approach the 
value a^^^Kp — arcosh [2 — cos(27rp/g)]. The condition 
r = for a transparent interface then gives the desired 
expression for the coupling constant, which after some 
algebraic manipulations can be written as 




The value 7^^^) = 7.8617*^^-' used in the numerical 
computations finally follows when Eq. (jB7|l is evaluated 
for 7^ ~ q'^'j^ and a lattice-constant reduction factor 
q = 5. 
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